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We report an updated measurement of the top quark mass obtained from pp collisions at \/s = 1.96 TeV at the 
Fermilab Tevatron using the CDF II detector. Our measurement uses a matrix element integration method to obtain 
a signal likelihood, with a neural network used to identify background events and a likelihood cut applied to reduce 
the effect of badly reconstructed events. We use a 2.7 fb — 1 sample and observe 422 events passing all of our cuts. We 
find mt = 172.2 ± 1.0 (stat.) ± 0.9 (JES) ± 1.0 (syst.) GeV/c 2 , or mt = 172.2 ± 1.7 (total) GeV/c 2 . 

1. INTRODUCTION 

The top quark is the heaviest known particle in the Standard Model. Its mass is an important parameter to be 
determined, both for its intrinsic interest, and because precision measurements of the top quark mass, in conjunction 
with the W boson mass, allow us to set constraints on the mass of the Higgs boson within the Standard Model. In 
this letter we describe a precision measurement of the top quark mass using a matrix element integration method. 
This measurement uses 2.7 fb _1 of data collected by the CDF II detector. 

We obtain a top mass measurement by integrating over unmeasured quantities in the matrix element using a 
quasi-Monte Carlo integration. This allows us to minimize assumptions made about the kinematics of an event, 
resulting in improved precision. The integration method yields a likelihood curve as a function of the top pole mass. 

The largest source of systematic uncertainty in our measurement is the jet energy scale (JES). To reduce our 
uncertainty due to this source, we introduce an additional parameter to our likelihood, Ajes 7 which allows us to use 
the information in the W decay to determine the JES. Ajes parameterizes the shift in JES in units of the systematic 
error for a given jet. Our likelihood is thus constructed as a 2D function of m t and Ajes! we then combine the 
likelihoods for all events and eliminate Ajes as a nuisance parameter to find a final top mass value. 

2. EVENT SELECTION 

At the Fermilab Tevatron, top quarks are predominantly produced in ti pairs, where the t decays into a W boson 
and a b quark ~ 100% of the time. The W can then decay into a charged lepton and a neutrino ("leptonic" decay) 
or a quark-antiquark pair ("hadronic" decay). We search for events in the "lepton + jets" channel, where one W 
decays hadronically and one leptonically. Thus, we analyze events with four high-energy jets (two from the b quarks 
and two from the hadronic W decay), at least one of which is required to be 6-tagged using a secondary vertex 
algorithm which identifies tracks displaced from the primary vertex; exactly one high energy electron or muon (from 
the leptonic W decay); and large missing transverse energy (from the neutrino). 

The principal backgrounds to our signal are events where a W boson is produced in conjunction with heavy flavor 
jets (bb, cc, or c), a W boson is produced with light jets which are mistagged as 6-jets, and QCD events not containing 
a W where the W signature is faked. Overall we expect 105.7 ± 42.0 background events in our observed 494 candidate 
events. Table Q] shows our expected backgrounds. 

3. MATRIX ELEMENT METHOD 

We calculate a two-dimensional likelihood as a function of m t and Ajes by integrating over the matrix element 
for ti production and decay over the unknown parton-level quantities, using transfer functions to connect these with 
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Table I: Expected backgrounds for the W+4 tight jet sample used. 



Background 


1 tag 


> 2 tags 


non-W QCD 


20.0 ± 17.3 


0.8 ± 1.6 


W+lignt mistag, diboson, or Z 


27.7 ± 5.6 


1.1 ± 0.1 


W+heavy (bb, cc, c) 
Single top 


45.0 ± 37.8 
3.8 ± 0.2 


6.0 ± 5.0 
1.2 ± 0.2 


Total background 
Predicted top signal 


96.5 ± 36.8 
259.4 ± 33.6 


9.2 ± 5.2 
98.8 ± 16.0 


Events observed 


389 


105 



the measured jets. Our overall likelihood formula is: 

1 1 24 

L(y\ m t ,Aj ES ) = — ( — t-tt t -^w i L i (y\ m f ,A JES ) (1) 

N(mt) A(m t , A JES ) ~— J 

with 

U(y | m t , A JES ) - J f{Zl) F p Z2) TF(y | f,A JES ) \M{m u x)\ 2 d$(x), (2) 

where x denotes the parton-level quantities, y denotes the quantities measured in our detector, M is the matrix 
element for ti production and decay, f(z) is the parton distribution functions (PDFs) for the momenta of the two 
incoming particles, FF is the flux factor normalizing the PDFs, N{m t ) is a normalization factor, A{mt, Ajes) is 
an acceptance factor to correct for the effect of the event selection criteria, and <& is the parton-level phase space 
integrated over. The integral is evaluated for each of the 24 possible jet-parton assignments and then summed 
with appropriate weights corresponding to the probability that a given jet-parton assignment corresponds with the 
observed 6-tags. We integrate over a total of 19 variables. In order to perform this integral in a practical amount 
of time, we employ quasi-Monte Carlo integration [l|, which uses quasi-random sequences. These sequences provide 
more uniform coverage of the phase space, resulting in faster integral convergence than with normal Monte Carlo 
techniques. 

We use a neural network to identify events likely to be background, and subtract out their contribution to the 
total likelihood by estimating the average contribution for background events from Monte Carlo. We also consider 
the effect of events which we call "bad signal" . These are events which contain an actual ti decay, but where the 
final observed objects in our detector do not come directly from ti decay (due to extra jets from initial or final state 
radiation, W — > r decay, or other causes). To reduce the effect of these poorly- modeled events, we apply a cut of 10 
to the peak of the log-likelihood curve. In Monte Carlo simulation, this cut eliminates 20% of "bad signal" events 
and 27% of background while retaining 97% of our good signal events. 

We test and calibrate our measurement using pythia Monte Carlo events over a variety of input m t and Aj E g 
values by performing pseudo-experiments. Using the results of the pseudo-experiments, we obtain a final set of 
calibration constants for our measured top mass and statistical uncertainty. Figure Q] shows the results of our Monte 
Carlo testing. Figure [2] shows the effect of the likelihood cut on Monte Carlo events. 

4. RESULT 

We have 494 events passing our initial selection cuts, of which 422 events pass the likelihood cut as well. With 
these 422 events, we measure: 

m t = 172.2 ± 1.0 (stat.) ± 0.9 (JES) ± 1.0 (syst.) GeV/c 2 = 172.2 ± 1.7 (total) GeV/c 2 
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Figure 1: Results using Monte Carlo events to test and calibrate our method. Left: measured mass vs. input mass. Right: 
pull width vs. input mass. 



Figure [2] shows the final contours of Tex, 2-cr, and 3-cr statistical uncertainty around the measured value. The total 
result attains a precision of better than 1% in wit- 

Our main sources of systematic uncertainty are from the Monte Carlo generator used for our calibration and 
testing (0.5 GeV/c 2 ), the residual JES uncertainty resulting from variation of the individual sources of our total 
JES uncertainty (0.5 GeV/c 2 ), uncertainty from the modeling of the jet energy scale for 6-jets (0.4 GeV/c 2 ), and 
uncertainty in the background model (0.4 GeV/c 2 ). We also have smaller uncertainties from initial-state and final- 
state radiation, lepton Pt measurement, pileup events, calibration, and PDFs (0.3, 0.2, 0.2, 0.1, and 0.1 GeV/c 2 , 
respectively), for a total of 1.0 GeV/c 2 . 
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Figure 2: Left: Value of the log-likelihood curve at its peak for good signal, bad signal, and background events in Monte Carlo. 
The dashed line shows the cut at 10 used. Right: Contours of the 2D likelihood distributions obtained with our final data 
sample. The contours shown correspond to a statistical uncertainty of 1, 2, and 3 a. 
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